Dissipation and supercurrent fluctuations in a diffusive NS ring 
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A mesoscopic hybrid Normal/Superconducting (NS) ring is characterized by a dense Andreev 
spectrum with a flux dependent minigap. To probe the dynamics of such a ring we measure its lin- 
ear response to a high frequency flux, in a wide frequency range, with a multimode superconducting 
resonator. We find that the current response contains, beside the well known dissipationless Joseph- 
son contribution, a large dissipative component. At high frequency compared to the minigap and 
low temperature we find that the dissipation is due to transitions across the minigap. In contrast, at 
lower frequency there is a range of temperature for which dissipation is caused predominantly by the 
relaxation of the Andreev states' population. This dissipative response, related via the fluctuation 
dissipation theorem to a non intuitive zero frequency thermal noise of supercurrent, is characterized 
by a phase dependence dominated by its second harmonic, as predicted long ago [T] [2] but never 
observed so far. 

PACS numbers: 



A phase coherent non superconducting conductor (N) 
connected to two superconductors (an SNS junction) 
gives rise to the formation of Andreev states (AS) 
which are coherent superpositions of electron and hole 
states confined in the N metal and carry the Josephson- 
supercurrent [3]. They strongly depend on the phase 
difference (p between the two superconductors. The 
quasi-continuous Andreev spectrum of a diffusive metal- 
lic wire exhibits a phase modulated induced gap, the 
minigap, which closes at odd multiples of 7r (for per- 
fectly transmitting contacts), and can be approximated 
by 2E g ((p) = 2E g {Q)\costp/2\ [H]. In the case of a 
wire longer than the superconducting coherence length, 
E g (0) ~ 3Eth, with Eth the Thouless energy H/td, and 
td the diffusion time across the N wire. Whereas most 
investigations of diffusive SNS junctions rely on non lin- 
ear transport measurements of current voltage curves, 
switching current and ac Josephson effect [7J [5] , only few 
experiments probe the Andreev spectrum at equilibrium 
in a phase biased configuration with NS rings threaded by 
an Aharonov Bohm flux. These include the tunnel spec- 
troscopy of the minigap [9] and the measurement of the 
flux dependent Josephson supercurrent using SQUID [TO] 
or Hall probe magnetometry [llj . Beyond probing the 
equilibrium AS spectrum in a static magnetic flux, the 
investigation of the dynamics associated to this spectrum, 
a far more complex question, has been addressed exper- 
imentally only recently 12J. Theoretically, it was pre- 
dicted that, in contrast to tunnel Josephson junctions 
[13] and because of the smallness of the induced gap, 
SNS junctions should exhibit low frequency supercurrent 
fluctuations at equilibrium [2]. According to the fluctu- 
ation dissipation theorem, in the linear response regime, 
such equilibrium fluctuations lead to a dissipative cur- 
rent under an ac flux excitation [T3][TS]. In this Letter, 
we present the linear current response of a phase biased 
NS ring and account for both the non dissipative compo- 
nent and the more surprising dissipative one, over a wide 



frequency range. We identify two fundamental processes 
leading to dissipation, the microwave-induced transitions 
across the minigap, and the energy relaxation of Andreev 
level populations. 

To this end, we couple a NS ring to a superconducting 
resonator, and phase bias it with a dc Aharonov Bohm 
flux and a small ac flux <5$ w at the resonator's eigenfre- 
quencies w. The linear current response 51^, is charac- 
terized by the complex susceptibility = 51^/8^^ = 
x'(uj)+ix"{u) = iujY(uj), where Y(w) is the NS ring's ad- 
mittance. This susceptibility is extracted from the varia- 
tions of the resonator's eigenmodes (frequency and qual- 
ity factor) . The ring's dissipationless response is deduced 
from the periodic flux variations of x', whereas the dissi- 
pation corresponds to \" . 

A first experiment [12] found a large dissipative re- 
sponse as well as a non dissipative one that differed no- 
tably from the adiabatic susceptibility, the simple flux 
derivative of the ring's Josephson current (the inverse ki- 
netic inductance). These results were partially explained 
by the theory of the proximity effect [15]. However, the 
shape of the flux dependences of \ did not vary with 
frequency (in the range explored), so that the different 
components of the ring's dynamical response could not 
be accessed. In particular, with the inelastic scattering 
rate much smaller than the lowest eigen-frequency, the 
dissipative response associated to the relaxation of An- 
dreev states could not be detected. 

In the present experiment, we report on a NS ring 
with enhanced temperature dependent inelastic scatter- 
ing rate 1 /t;„ (T) , thanks to a thin Pd layer at the NS in- 
terface, between the normal gold mesoscopic wire and the 
superconducting niobium loop. The higher inelastic scat- 
tering rate, combined to a broader frequency and temper- 
ature range, lead to the identification of the two funda- 
mental contributions to the supercurrent relaxation. At 
frequencies above the inelastic scattering rate, dissipation 
is due to microwave-induced excitations across the mini- 
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FIG. 1: a. Inset : Scanning electron micrograph of 
the NS ring fabricated to explore the dynamics of Andreev 
States. a.NS ring is connected to a Nb multimode resonator 
(meander in the top) by W wires deposited by a FIB. b. 
Flux dependence of \' a t several resonator eigenfrequencies 
(/o = 190 MHz, fi = 560 MHz, f s =2 GHz) and T = 1.2 K. 
At fo, x' is barely distinguishable from a cosine (solid line). 
Whereas the amplitude 8x^-0 = x'(^)~ x'(0) is frequency in- 
dependent, its temperature dependence is that of the switch- 
ing current I s . Note the appearance of a local maximum 
around tp = at high frequency, c. Temperature dependence 
of the switching current of the control sample I a (circles) and 
of |J j^-Sx'tt-o measured in the NS ring (squares) with their 
fits according to ref.[8] (respectively dashed and solid lines). 
Best fit yields Eth = 71 mK±5mK corresponding to 1.5 GHz 
for the NS ring. 



gap. In the opposite regime of lower frequency, which 
could not be reached previously, dissipation is due to the 
relaxation of Andreev level populations. 

Here we reveal this second contribution, proportional 
to the sum of the squared Andreev level currents. Ac- 
cordingly we measure a response whose period in flux is 
nearly 7r periodic. The extra cusps we find at odd multi- 
ples of 7T reflects the closing of the minigap. This charac- 
teristic phase dependence, which is precisely that of the 
low frequency, thermal supercurrent noise, is in complete 
agreement with theoretical predictions formulated long 

ago m mug. 

The experimental set-up is shown in Fig. la, the res- 
onator consists in a double meander line etched out of a 1 
micron thick niobium film sputtered onto a sapphire sub- 
strate. The NS ring connects the two lines at one end of 
the resonator, turning it into a A/4 line with a fundamen- 
tal frequency of 190 MHz, and harmonics 380 MHz apart. 
A weak capacitive coupling to the microwave generator 
preserves the high quality factor of the resonances, which 
can reach 5 10 4 up to 10 GHz. The NS ring is fabricated 
by electron beam lithography. The Au wire (4 micron 
long, 0.3 micron wide and 50 nm thick) is first deposited 



by e-beam deposition of high purity gold. The S part 
is deposited in a second alignment step by sputtering of 
a Pd/Nb bilayer (6 nm Pd, 100 nm Nb). The resulting 
uncovered length of the Au wire is 1/im. The ring is con- 
nected to the Nb resonator in a subsequent step, using 
ion-beam assisted deposition of a tungsten wire in a fo- 
cused ion beam (FIB) microscope. This process creates a 
good superconducting contact between the resonator and 
the Pd/Nb part of the ring. The 6 nm-thick Pd buffer 
layer ensures a good transparency at the NS interface, 
as demonstrated by the amplitude of the critical cur- 
rent measured with dc transport measurements on con- 
trol SNS junctions fabricated simultaneously (Fig.lc). It 
also enhances the inelastic scattering rate because of Pd's 
spin- wave like excitations (paramagnons) [TBHIH]- Con- 
sidering that the phase coherence time extracted from 
weak localization measurements on a 6nm thick Pd thin 
film [16], was of the order of 0.3 ± 0.1 ns at 1 K, which 
is longer than the estimated diffusion time td= 0.1ns 
through the Au wire between the S contacts, we do not 
expect a reduction of the critical current as confirmed by 
measurements in the control samples. 

The quantities we measure are the variations with 
dc flux of the resonator's quality factor and eigen- 
frequencies 8Q($) and <$/($). They are simply related 
to the oscillating phase dependent part of the complex 
susceptibility, characterized by xf(<p) and x"^)) where 
the superconducting phase tp is related to the flux thread- 
ing the ring by ip — — 27r<I>/$o where <&o = h/2e is the 
superconducting flux quantum. The relation reads |12j : 
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The coupling inductance Lc — 9 ± 2pH is due to the 
S part of the N/S ring; C —0.3 [iH is the inductance 
of the resonator. These expressions are valid at tem- 
peratures such that the kinetic inductance of the SNS 
junction is larger than the ring's geometrical inductance 
(outside this range screening of the applied flux, both dc 
and ac, needs to be considered [12]). This sets the lower 
limit to the temperature, so that experiments were con- 
ducted between 0.4 and 1.5 K. The frequencies probed 
ranged between 190 MHz and 3 GHz. 

We find spectacular variations of both the amplitude 
and shape of x' an d x" as frequency and temperature are 
changed. At the lowest frequencies and highest tempera- 
tures investigated (see Fig.l and Fig. 2a) the dissipation- 
less x'if) i s weu described by a pure 27r-periodic cosine, 
as expected for the adiabatic susceptibility xj = 91 j /d& 
of the Josephson current which is purely sinusoidal at 
these moderately high temperatures, much larger than 
Eth [S HI]- As shown in Fig.lb,lc, The amplitude 
&x'tt-o — x'M — x'(0) perfectly reflects the expected, 
roughly exponential, decay of the Josephson critical cur- 
rent Ij(T) = Ij(0)exp(-k B T/3.6E Th ) [H], that was also 
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FIG. 2: Evolution of x' (top) and x" (bottom) phase depen- 
dence with temperature at fi = 560 MHz < Erh/h (left) and 
f B = 1GHz > E Th /h (right). X '{f = t) and x "(v = t) in- 
crease with decreasing temperature. x'(.<P) an< A X" (v 3 ) strongly 
differ at low frequency and high temperature whereas they 
are similar with a shape reminiscent of the minigap at high 
frequency. The solid line in d is proportional to the mini- 
gap dependence \cos(tp/2)\. Curves have been shifted so that 
X'(<p = tt) - X '(<P = 0) = and x"(v = 0) = 



measured in the control wire. We find for both samples 
Exh = 71 ± 5mK which corresponds to to — 0.1ns. 
In contrast, the dissipation, characterized by x"^)) is 
nearly tt periodic (see Fig. 2b) at the largest temperatures 
investigated and acquires a strong 2ir periodic component 
at lower temperature. When increasing the frequency, 
(Fig. lb and Fig. 2a and 2c ) x'if) contains additional 
harmonics, with peaks at odd multiples of tt and more- 
over a local maximum at (p = 0, mod [27r] for the highest 
temperatures. On the other hand, at 2 GHz (Fig2.c and 
2.d ) and low temperature, x'(v) an d x"(v) have iden- 
tical shapes, with peaks at tt, mod [2tt], reflecting the 
underlying minigap that varies like \cos(tp/2)\. 

In the following we exploit this complex evolution of 
x(<p) with frequency and temperature to extract the dif- 
ferent mechanisms at work in the dynamics of Andreev 
states. To this end we make use of theoretical predictions 
[15j based on Usadel equations, and recent numerical sim- 
ulations [T!5] inspired by the analysis of the ac response of 
normal mesocopic rings [201422J . The response function 
of a NS ring has been shown to contain three contribu- 
tions: x — XJ + Xd + Xnd- The adiabatic, zero fre- 
quency, Joscphson contribution x,J is purely real and is 
the derivative of the Josephson current dlj / <9$ . The sec- 
ond contribution, the diagonal susceptibility xd, is the 
first non adiabatic, frequency dependent contribution. It 
describes the Debye-like relaxation of the (phase depen- 
dent) thermal populations f n {<P) of the Andreev states, 
with a typical inelastic relaxation time Tj„ according to 



the simple model proposed for the dynamics of persistent 
currents in normal rings 20Tl22j : 
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where the square of i n = — <9e n /<9<E>, the current car- 
ried by the n-th Andreev level of energy e„($), ap- 
pears. Finally, the non-diagonal contribution xnd de- 
scribes quasi-resonant microwave-induced transitions be- 
tween two Andreev levels, involving (in contrast with 
Xd) non diagonal matrix elements of the current opera- 
tor [TH]. The contribution x'nd to tne phase dependent 
susceptibility dominates when oj > 1/to 3> 1/Tj n , as 
in Fig. 2d. x' an d x" then have similar shapes which 
follow approximately the minigap with peaks at tt and 
a I cos (ip/2) \ dependence [19]. This high frequency regime 
was the only one accessed in the previous experiments on 
Au wires directly connected to W superconducting wires. 
In those experiments the energy relaxation time, limited 
by electron electron interactions, of the order of O.lyus, 
was very long due to the superconducing contacts [T2| [23], 
Therefore those measurements were always in the regime 
WT m ^ 1 where x'r> 1S negligible, (Eqj2]). In contrast, the 
Pd layer beneath the Nb contacts in the present sam- 
ples considerably reduces the inelastic scattering time, 
leading to a substantial contribution of x'd f° r the res- 
onator's first five eigenfrequencies. We now focus on this 
contribution analyzed in Fig. 3 and Fig. 4. 

We first present the predicted flux dependence of xd, 
given by the function F 
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which reads in the continuous spectrum limit _F($,T) = 
JdeJ|($,e)/ [4fc B Tp(e)cosh 2 (e/2fc i3 r)]. Here J s and 
p are respectively the spectral current and the density 
of states of the SNS junction. This function was intro- 
duced by Lempitsky [T] to describe the I(V) character- 
istics of SNS junctions, and was calculated numerically 
using Usadel equations by Virtanen et al4l5 . At large 
temperature compared to Exh, F(<f>,T) can be approx- 
imated by the following analytical form: Fu(ip,T) cx 
[(-tt + (tt + v)[27r])] sin(y>) - I sin(<^)| sin 2 (^/2)/7r. It is 
dominated by its second harmonics with in addition a 
sharp linear singularity at odd multiples of tt (see Fig. 4). 
This is due to the dominant contribution of Andreev lev- 
els close to the minigap whose flux dependence is singular 
like in a highly transmitting superconducting single chan- 
nel point contact [2]. We now show that the particular 
flux dependence of the function Ftj can explain the exper- 
imental data of Fig.l and 2. We first follow the frequency 
dependence of the amplitude of x'd^) = x'iv) ~ Xj(v) 
at fixed temperature, and check that the shape of Xz?(</0 
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FIG. 3: a. Experimental phase dependence of x'd a t 
fi = 560 MHz and (from top to bottom) T = 580 mA", 
650 mK, 720 mK, 820 mK, 1000 mA", 1200 mA". In agreement 
with theory [15], the position of the maximum slightly shifts 
to lower values with decreasing temperature, b. Frequency 
dependence of &Xd '■ the maximum of — x'oi'p), a t different 
temperatures (triangles) compared to the theoretical predic- 
tion from Eq.3. c. Temperature dependence of extracted n n 
(see text) . Solid line is a T~ 3 law. Inset: frequency depen- 
dence of 5x'r> & t 1.15K (circles) compared to the theoretical 
prediction from Eq.3. 



does not change with frequency and is the same as that 
of x'b, as predicted for the temperature and frequency 
regime where the contribution of x'nd can be neglected. 
As shown on Fig. 3b it is then possible to fit the frequency 
dependence of the amplitude of x'd if) by the expected 
(wTj n ) 2 /(l + (ujT in ) 2 ) and determine the characteristic 
time Ti n for several temperatures according to Eqj2] . We 
find values of r% n varying between 0.2 and 1 ns, quite 
similar to what was deduced from weak localisation mea- 
surements in Pd thin films [16] . Moreover the power law 
decrease of r in (T) in T~ 3 , is in reasonable agreement 
with what is expected for paramagnons which constitute 
the dominant inelastic scattering at low temperature in 
Pd which is close to a ferromagnetic transition. It is also 
interesting to note that our results can be described by a 
single inelastic time, independent of ip, whereas a phase 
dependent r in is expected for electron phonon collisions 
in SNS junctions [23]. This is probably due to the fact 
that temperature is larger than E g (0) in our case. 

A similar analysis can be done on x" , the quality of the 
calibration is however not as good as on x' ■ Moreover we 
still lack a good analytical prediction for xnd i^P) > which 
gives a large contribution to x"(<^) & t l° w temperature 
and high frequency. We have overcome this difficulty by 
subtracting for frequencies larger than 1.7 GHz the flux 



FIG. 4: Normalized experimental phase dependences at 
T — 1.2 K (open symbols) compared to Fu (solid). Data are 
X" (squares), \x'd\ — XJ ~ x! (circles) at /i = 560 MHz and 
X" (triangles), \x'd\ — XJ ~ X (diamonds) at — 1.35 GHz. 
At high temperature and for / < Erh/h, experimental XJ~ x' 
and x" are found to be in very good agreement with theoret- 
ical prediction of Usadel equations Fu ■ 



dependence of x'nd estimated from the high frequency 
data (2.8 GHz). The resulting amplitude 6x'b(u) agrees 
with the expected frequency dependence in wt;„/(1 + 
(cjTj„) 2 ) as shown in the inset of Fig. 3c. One can also 
compare the independently measured flux dependences 
°f x' ~ Xj an d x" with theoretical predictions from the 
Usadel equations, Fjj(§,T). This is done in Fig. 4 for 
several frequencies and a good agreement is found. 

With this set of experiments, we have thus shown that 
the frequency and temperature dependences of the re- 
sponse function of NS rings in a time dependent flux are 
consistent with a simple Debye relaxation model of the 
population of the Andreev levels. Using fluctuation dis- 
sipation theorem one can estimate the related thermody- 
namic current noise as: 



2 k B T X ' D {u) 2 
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1 + (WT,„) 

(4) 

The measurement of the ac current linear response of 
an NS ring to an ac flux thus reveals two fundamen- 
tal mechanisms contributing to dissipation at finite fre- 
quency. One of them, predominant at high frequency 
and low temperature, describes the physics of microwave 
induced transitions above the minigap. We have clearly 
identified and characterized the second cause of dissipa- 
tion, the thermal relaxation of the populations of the 
Andreev states. It is described by an inelastic rate which 
is extremely sensitive to the nature of the NS interface. 
This dissipative response is directly related to the low 
frequency thermal noise of the Josephson current, with 
a flux dependence proportional to the average square of 
the spectral (or single level) current, and can be pre- 
cisely described by theoretical predictions. These results 
show that linear ac measurements in a wide range of frc- 
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quency, close to equilibrium, reveal physical properties of 
SNS junctions that are not accessible by standard trans- 
port measurements dominated by non linear effects. The 
type of experiments presented here is uniquely suited to 
investigate more exotic systems, for instance with the 
normal diffusive wire replaced by a ballistic wire, leading 
to a discrete Andreev spectrum known to be extremely 
sensitive to spin orbit interactions. 
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